ON THE DISTRIBUTION OF THE INTEGRAL OF THE 
EXPONENTIAL BROWNIAN MOTION 



LEONID TOLMATZ 

Abstract. The density distribution function of the integral of the exponential 
Brownian motion is determined explicitly in the form of a rapidly convergent 
series. 



1. Introduction 

Let B{s), < s < t denote a standard Brownian motion. 
In the present paper we consider the functional 

fe'^<^^^ds. (1) 
Jo 

This functional naturally appears in mathematical finance in connection with 
Asian options and in various other applications. Understanding of a functional is 
stipulated by computability of its distribution, and this problem has attracted at- 
tention of researchers since 1992 when its study was initiated by M.Yor. For details 
and further references see Matsumoto and Yor (2005a), (2005b), Dufresne (2000), 
Linetsky (2004), Schroder (2006). 

In the present article we obtain a new representation of the density distribution 
function f{X,t) of this functionl in the form of a rapidly convergent series. This 
series is given in Theorem 11 Section 8. It allows computations with any desired 
accuracy in the whole domain, except small values of the argument, but for these 
asymptotic expansions are available; these will be discussed elsewhere. 
These results were obtained by the author between 1997 and 1999, but remained 
unpublished. In 2003 they were announced on the author's web page with some 
additional results and were communicated at the AARMS workshop on Financial 
Mathematics in St. Johns in 2003. 

The method applied readily leads to similar series for the cumulative distribution 
function and pricing of Asian options. Our approach is based on the classical 
method of M.Kac. 

Notation and Conventions. 

C(^) is Riemann's Zeta function, 

iFi(a;7; z) is the confluent hypergeometric function. 
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Hm{x) arc Hcrmite polynomials, 
C = 0.57721 ... is Euler's constant, 
[x\ is the integer part of x, 

F{X, t) and /(A, t) are the distribution and density distribution functions of func- 
tional (1). 

2. Double Laplace Transforms and the Corresponding Differential 

Equation. 

Let 



F{X,t) := p| ^* e^^Wds < A I 



(2) 



be the distribution function of functional (1) and let u{p, q) denote its double 
Laplace-Stieltjes transform: 

/■+00 r+oo 

u{p,q)= / e-P^-i'dxF{X,t)dt. (3) 

Jo Jo 

A direct application of the classical result of Kac (1949), see also Rosenblatt 
(1951), shows that the function u{p, q) can be recovered in the form 

/ + 00 
u{p, q, x) dx, (4) 
-oo 

where u{p,q,x) is a solution to the following second order ordinary differential 
equation 

u"{p, q, x) — 2{q + pe^^) u{p, q, x) = 0, where x ^ 0, (5) 
for which we seek a continuous solution subject to conditions 

u{p, q, x) ^ 0, as a; — > ±oo, (6) 

\u'{p,q,x)\< M, x^O, (7) 

u'{p,q,0+)-u'{p,q,0-) = -2. (8) 

The distribution F{X,t) and its density f{X,t) := Fx{X,t) can be recovered from 
u{p, q) by double Laplace inversions: 

^(^'^) = 7^r-W / '^'i^'' -u{p,q)eP^dp (9) 

/(^'*) = 7;r-W / '^'l^' l u{p,q)eP^dp (10) 
with any a > 0, 6 > 0. 

Theorem 1. The general solution of the equation 

y"{x)-2{q + pe^'')y{x)=0 (11) 

^5 given by 

y{x) = ci/^(V2^e^) + C2K^{^/2^e-), (12) 
where Iv{x) and K^{x) are modified Bessel functions of the first and second kind. 
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Proof. Equation (5) is a particular case of the equation 

x'^y" - 2axy' + [(b^e^^^ - v'^)c^x'' + a{a - l)]y = 0, 

where a = 0, c = 1, 6^ = -2p, v"^ = 2q, see Kamke [1959] 2.163 (23), and its 
solution is 



y = x-^Z,{be'^) = Z^(v^e^) = Z^^{±i^e-) = 

□ 

Theorem 2. The solution of equation (5) subject to conditions (6) - (8) is of the 
form 

uiv a X) = / 2/^(V2^)if^(V2^e^), z/x > 0, 
' \ 2K^^{^)I^^{^e%ifx<^. 

Proof. According to Theorem 1 we seek a solution of the form 

ci7^(^e-) + C2K^^{^e-), if a; > 0, 
c'J^qiV^e-) + 4/^^(^/2^e-), if a; < 0, 

where constants ci, . . . , to be determined by conditions (6) - (8). 

Since K^{x) — » oo as a; — » 0+ and it decays as x ^ oo, while I^ix) — » as a; — + 0+ 

and Iv{x) — > oo as a; — > oo, in order to satisfy condition (6) we set 

ci = 4 = 0, 

hence 

\ c'i/v2^(%/2^e-), ifa;<0, 
and this function obviously satisfies (7). 



u{p, q, x) 



It remains to show that the choice 

ci = 27^(^), 

generates a continuous solution that satisfies (8). 
The continuity follows from 

u{p, q, 0+) = u{p, q, 0-) = 2 V2p)^V2^(\/2p). 
By a straightforward computation we obtain: 
u'{p,q,0+)-u'{p,q,0-) = 

We observe that the expression in the brackets is up to the sign the Wronskian: 

-W{K4z),L{z)} = --, 

z 

see Olver (1974) p. 251 (8.07), where z = i/2p, = -/2g, and this yields 

u'{p,q,0+)-u'{p,q,0-) = -2. 

□ 
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In order to recover the density distribution /(A, t) from u{p, q, x) we first invert the 
latter with respect top (Section 3), then integrate the result on x (Section 4 ), and 
then invert with respect to q in the following sections. These steps involve inter- 
change of integrations which is easily justifiable by simple estimates and Fubini's 
theorem, see for example Tolmatz (2000), (2003). 



3. Laplace Inversion with Respect to p 
The following theorem provides a Laplace inversion of u{p, q, x). 
Theorem 3. 

MA,g,.) ■■-^J^_^^ uip,q,x)e^ *=A^"H"^^J'^Uj-(,4) 

Proof. It suffices to verify that the Laplace transform of the right-hand side of (14) 
results in u{p,q,x), that is 



J w{\, q, x)e~P^ '^^^ j A V 2A~ j ^^\\) ^"^^ '^^ ^ ^' 

To this end we make use of a table of Laplace transforms, see for example Ditkin 
and Prudnikov (1961) 9.351: 

r^'P ( " ^) (^) '^^ = ^^-^ + ~ ^^^16) 

where > 3?/3 > 0. 
We apply (16) with 

r a + /? = + 1, 
1 a- (} = 2e^, 



that is 



and obtain 



a = i(e^ + l)2, 
/3=i(e--l)^ 



^1 \/2e^, if a; < 0. 



(18) 
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Substitution of (17)-(18) in (16) yields 



2A 



-p^ dX 



\ 27^(V2^)K^(V2^e-), if a; > 0, 
\ 2K^{V2^)I^{^e-), iix<0. 
u{p,q,x). 



□ 



Remark. Alternatively, we can prove the above theorem by a straightforward deriva- 
tion without "guessing" the function in (14), however that proof involves tedious 

computations. 

Remark. w{X, q, x) is a Laplace transform of the joint density distribution of random 
variables J* e^^^'Us and B{t). 

4. The Integral Jwdx 
In this section we compute the integral 

/•oo 
u{X,q,x) dx. 
-oo 

The following proposition will be used in the proof of the next theorem. 
Proposition 4.1. Let 3?a > 0, 5?(/i + > -1, > 0. Then 

e-"""^ I^{l3x) dx = 



J 

Jo 



(19) 



/3''r(i^. + iM + i) 



2''+iQ!5(''+''+i)r(i/-i-i) 

Proof. By Gradshtein and Ryzhik (1963) 6.631 1. 

e-"""^ JAj3x) dx = 



; v + I; — 

2 ' '4a 



POO 

Jo 

2^+iaiii^+^+i)r{i' + 1) ^ 



(20) 



Fi 



1 

4a 



The integral in (20) is a function of the real parameter that admits an analytic 
continuation into the domain D := {(3\ arg/3 < tt} in the complex /3-planc. 
Indeed, for any fixed x > J^iPx) is analytic in D, while the asymptotics of Jv{z) 
in the domain | arg^j < tt according to Gradshtein and Ryzhik (1963) 8.451 1. is 



and this implies 



2 

TTZ 

2 

TTZ 



COS 



%-J)(i + o(i/|.n) 



as \z\ 



cos(/3x-|z.-|)(l + 0(l/|/3n) 



as 



oo, 



which ensures that for any (3 & D the integral in (20) converges and is differentiable 

on (3. 

The function in the right-hand side of (20) is also analytic in D due to the analycity 
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of iFi(a;7;z) in the whole z-plane and since fi" is analytic in D. Therefore (20) 
holds also on the imaginary /J- axis, and we get 

Jo 2-+iai(''+''+i)r(i/ + l) V 2 ' ' 4a y 

The proof is completed by substituting 

in the previous relationship and obvious simplifications. 



(21) 



□ 

Theorem 4. 

Proof. In the integral 

wiX,q) = J u{X,q,x)dx = J ^ ( 2A~ ) T J 

we perform the change of variables 

y = e"^, 
X = Iny, 

dx — —, 

y 

and obtain by Proposition 4.1 with fj, = —1, a = 1/(2A), (3 = 1/X, u = ^J2q 

dx = (22) 



»(A,,)4e-"(-.£exp(-^)/^( 
le-./,«,£ie-.V™,^(i),,= 



r(^) ^(V^. , 1 



□ 



In order to compute the distribution density /(A, t) of the functional Jq e^^^^^ds 
we seek to invert w(A, g), that is to compute the integral 

-I /'C+ioo 

fM = ir- w{\,q)e'i'dq= (23) 

ZTT? Jc—ioa 

These computations will be performed in Sections 5-8. 
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5. A Series for w{X, q) 



Theorem 5. 



where 



»(A,<,) = i;e-"™f;l.,(A,<,)pl, 



n=0 



Wn{X,q) := (2A)-v^/2 



r(^ + 



r(^ + n+l)' 
Proof. The proof follows from Theorem 4, the expansion 

iFi a; 7; z = 1 + - - + ' 

7 1! 7(7 + 1) 2! 

where a = y/2q/2, 7 = i/2g +1,2= 1/(2A) and the identity 

r(a + n) = r(a) a ... (a + n - 1), 

where a = ^/Sg/S, \/25 + 1. 



(24) 



(25) 



□ 



In the remaining part of this section we shall justiiy a term-wise inversion of the 
series in (24). 

The following elementary lemma is given without proof. 



Lemma 5.1. Let \ a.rgz\ < Tr/2. Then for any s > 

z/2 + s 



z + s 



< 1. 



Corollary 1. For \ a.rgq\ < n, n = 0, 1, . . . holds 

/2g/2 + n 



2q + n 



Lemma 5.2. Let \^q\ > 0, n = 0, 1, . . . 

Then 

r(v^/2 + n 



r(V2^ + n+l) 



< 



2q\ 



r(v^/2) 



r(V2g) 



(26) 



(27) 



(28) 



Proof The proof follows by elementary transforms and Corollary 1: 



r(v^/2 + n) 



r(^ + n+l) 



1 



(72^/2 + n - 1) . . . iy^/2) r(V2^/2) 



< 



< 



r(v^/2) 



^ + n r(v^) 



< 



2q\ 



r(V2g) 



Lemma 5.3. Let \^q\ > 0, n = 0, 1, . . . . 



r(v^/2 + n) 



r(^ + n+l) 



< 2v^. 



1 



^1 



I2-V29 I 



r(^/2 + 1/2) 



□ 



(29) 
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Proof. The proof follows by the previous lemma and Legendre's formula 

r(2a) = ^22"-ir(a)r(a + 1/2), (30) 



which we apply in the form 



-2V^2-'^—^—-, (31) 



r(2a) " r(a + l/2)' 

where a = ■^/2g/2. □ 

Lemma 5.4. Let ^(t) = c + ir, where c > and — oo < r < oo is a parameter. 

Then 



r(V2^/2 + n) 



<C'|^|-«\/M^ (32) 



r(V2g(7)+n + l) 
w/iere C, a > are some constants not depending on n. 

Proof. The proof readily follows from the previous lemma and Stirling's formula 

r{z)~e-'z'l^\ (33) 

as 2; — > oo in the sector | arg2;| <'it — 6, which implies 

1 



= 0(|r|-"V^) (34) 



r(V2(c + ir)/2 + l/2) 
as r — > ±00. 

□ 

The following theorem shows that the series in Theorem 5 can be inverted term by 
term. 



Theorem 6. 



where 



/(A,t) = ^e-V(-)X:-/4A,t)^, (35) 



^ pc-\-ioo 

U\t):=— WniX,q)e''*dq, (36) 

ZTTi Jc—ioo 



where 



anrf i/iis series converges absolutely for any fixed X,t > 0. 
Proof. The previous lemma implies absolute summability of 



n!^ r(V2^ + n + l) (2A)" 
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as a function of variables t, n with any fixed A, t > in the sense 



dr < OO, 



(2A)" 



(39) 



r(v/2g(^ + n + i) 

and the proof follows by Fubini's theorem. □ 

6. Transforms of ^^(A,^) 
We wish to prepare the integrals 

-| pc+ioo 

U\t) = — Wn{\q)e'^'dq (40) 

ZTTi J c—ioo 

for an explicit computation. 
Proposition 6.1. 

1 f-c+ioo 

fn{\t) = — Wn{X,q)e'^'dq= (41) 

Zm 



e 



C — IOC' 

,t/2 ^ fC+ioo 



2A 27rz 



Jc-ioo r(2Vi + n + 2) 



Proof. The proof follows by the reparametrization = +1, q = 2z + 2\fz + 
1/2, d9= (2^+l)/(^)rf2. □ 

We shall need auxiliary functions 

^„(^) ^^2^^r(^ + »+ 1/2) 2^ + 1^ (42) 
^"^ ^ r(2u; + n + 2) u; ' ^ ^ 

where n = 0, 1, 

By making use of the properties of the Gamma function as in the previous section 
we can write 

while for n = 1, 2, . . . 

^^(^)^2^.r(- + " + 1/2)2^ + 1^ 

^"^ ^ r(2w; + n + 2) w ^ ' 

^2-(»+i) ^ (w; + n-l/2)-...-(«; + l/2) 2u; + l 1 



«; + n/2 + l/2 (w + n/2) • . . . • (w + 1/2) w r(«; + 1) 

^/^2-^"+^^ ^ / " .. + j-l/2 \ 2z. + l 1 

^ «; + n/2 + 1/2 V^^^^i w + i/2 y w r(w; + l)' 

Consider the rational factor in (44) 

r iw) ■= ^AF2-("+i) 1 ( fr w + j-l/2 \ 2w + l 

Observe that due to cancellations of certain factors in the nominator and denomi- 
nator ro(w) is regular except the simple pole at w = 0, and r„(w), n > 1 is regular 
except the simple poles at w = 0, — 1, . . . , — [^^^] — 1- 

We wish to decompose r„(t«) into a sum of elementary fractions. 
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Theorem 7. 



[(n-l)/2] + l („) 

^"H= E ^> (46) 

fe=0 



(„) _ r(n + l/2) 



r(n + 2) ' 



(fc-1)! r(n-2fc + 2)' (^^^ 

w/iere fc = 1, . . . , [(n - l)/2] + 1. 

Proof. Since all the poles of r„(«;) are simple, the existence of decomposition (46) 

is a well-known fact. 
The relationship 

w„(w)r(w + 1) = 2''"'-^7 = (49) 

^"^ ^ ^ r(2u; + n + 2) w ^ ' 

^^.-(n+i) 1 (^ + n-l/2).....(^. + l/2) 2^. + ! ^ [("-^^1+^ ^ 

w + n/2 + 1/2 (w + n/2) •...•(«; + 1/2) w w + k 

shows that the coefficients aj."' are the residues of (^„(w)r(w+l) at w = 0, — 1, . . . , — [(n— 
l)/2] - 1, hence 

4") = Res^=.k{Mw)T{w + 1)} = i?es^^_ J22-I^;j±^±i/^?!^r(t« + 1)| = 

i {Iw + n + 1) w ) 

C r(n+l/2) :f^_n 
_ ) r(n+2) ' II « — U, 

\ ^n{~k)Res^=^knw + 1), if fc = 1, 2, . . . , [ii^^] + 1, 

f r(n+l/2) jfL_n 
^ I r(n+2) ' II «^ - 

1 ^„(-fc)((^, if fc = l,2,... + 
r r("+i/2) ;ft,-o 

^ I r(n+2) ' II - 'Jj 

1 9-2fc (-l)''^'' r(«-fc+l/2) •fjL_-| r, ["-11 II 

L ^ (fe-1)! r(n-2fe+2) ' ~ ^' • • • ' I 2 J 

□ 

For typographical reasons we introduce the notation 

a := 8Ae-2*. (50) 

Theorem 8. 

/(A,f) = ^— e-V(2A)e*/2y 1 /„(A,t)-i^, (51) 
^ ^ 2\V2X ;^n! -""^ ^ (2A)" ' ^ ^ 

w/iere 

/o(A, i) = / " rf^, (52) 

271-1 ic-ioo \/^r(v^ + i) 
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and for n = 1,2, 



[(n-l)/2]+l 



fn{\t)= a 



(n) 



fc=0 



2TTi 



C+ZOO 



(yi + fc)r(yi+i) 



dz, 



(53) 



where coefficients a^"^ are defined in Theorem 7. 
Proof. Only the case n > 1 needs a proof. By Theorem 7 we can write 

r(v/i + n+l/2) 2Vi+l _ _ ^_2^ r„(^/i) _ 

r(2yi + n + 2) ~ ^nWz)-^ r(vi+i)~ 

[(«-l)/2] + l („) 

fro ^^ + ^' 

substituting of which in (41) yields the proof. 

7. Auxiliary Results 
The following expansion is given in Gradshtein and Ryzhik (1963) 8.321 2: 



□ 



^ OO 

+1) /-^ ' 

' k=0 



r{z + i) 

where 

[ dn+l = :^J2'k=o(~^)''^k+ldn-k, 

whore ,si = C, s„ = C,{n) for n > 2. 
Theorem 9. 

|d„| =0(n-(i-^)"). 
Proof. By Cauchy's formula for Taylor coefficients we can write 



(54) 



(55) 



(56) 



dr. 



1 



1 



27ri Ji^i^i 2"+ir(^ + l) 



dz, 



and then estimate this integral by integrating over the circle z = rnB^'*', where 
= n/lnn : 



\dn\ 



1 

1 

2^ 



2"+ir(^ + i) 
2- # 



1 r' 



rn+^e'("+i)¥'r(r„e'V') 
/2 # 



1 

2^ 

< 



1 



2 (-Ztt 



27rr 



n+l 



d(/3 



|r(r„e*'^)| 



+ 



1 /■3/27r 
^X/2 



|r(r„e'v)| 27rrJJ+W./2 |r(r„e'v)|- 

Estimates for the latter two integrals can be obtained by estimating their integrands 
by Stirling's formula (33) with making use of Euler's formula 



zT{z)T{-z) 



SmTTZ 



(57) 
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in the second integral: 

J-TVI 



.11 |r(r„e''^)| 



= / O \^f7n 

f 



/2 |r(r„e^v)| 

7r/2 



e'-f sin (7rr„ 

J-7r/2 V 



and therefore 



o(l+7r)r„ r„ 

' n ]i 



0( e(l+'^)»/lnn( _!L 

Inn 



-n(l-l/lnn-l/(2n)) 



0(n-(i-^)"). 



□ 



Remark. By applying the saddle point method we can obtain an exact asymptotics 
for dn and thus improve the above estimates. 



Theorem 10. 



where 



8. Inverse Laplace Transforms 



^ ^ 2AV2A fb^n! -""^ ^ (2A)" ' 



n=0 



dz, 



(58) 



(59) 



anrf for n = 1, 2, . . . 



[(n-l)/2] + l 



fn{\t)= ^ a 



(n) 



fe=0 



(0+) 



27r«7-oo (v^ + A;) r(Vi + 1) 

(n) 



dz, (60) 



where a = 8Ae *, and the coefficients a). ' are defined in Theorem 7. 

Proof. The proof follows from Theorem 8 by pushing the integration contours to 
the left, which is justified by Cauchy's theorem and obvious estimates by Stirling's 
formula (33) with making use of Euler's formula (57). □ 

Lemma 8.1. 

27ri 



- / (yj)"-ia-^e2^*dz =o\ 



(61) 
(62) 
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where n oo, and a,t > and A: = 1, 2, . . . are fixed. 



Proof. Estimate (61) follows by integrating over the lower and the upper cut arg z = 
=F7r, an obvious reparametrization, elementary estimates and Stirling's formula: 



2tt Jo 27r Jq 



dx = 



n/2 



Estimate (62) follows in a similar fashion. 
Proposition 8.2. 

oo /— /•(0+) 

/o(A,i)= E^™^ / {yTzT-' a-^^ e^^' dz, 

m=0 



□ 



(63) 



and for n = 1, 2, . . . 



[(n-l)/2] + l .00 

/„(A,i)= ^ 4")^d„_L/ 



(04 



fc=0 m=0 

where all the series are absolutely convergent. 

Proof. The proof follows by substituting the expansion 



^/z + k 



a-^e^'^Uz, (64) 



oo 



r(Vi + i) 



(65) 



fc=0 



in the integrals (59) and (60) and term by term integration. 

The term by term integration is justifiable by absolute convergence and Fubini's 
theorem. The absolute convergence follows from estimates (56), (61) and (62): 



- 1 /■(0+) 

^/^^i-oo 



sfz + k 

Eo(„-<.-.)".,o{(|)"']< 



m=0 

oo 



/2 



m=0 



1 /•(°+) iVzY 



2m 



y/z + k 



a-^e^'Uz 



oo, 



where fc = 0, 1, . . . and < e < 1/2. 

Proposition 8.3. Let a he a real parameter and m = 0, 1, . . 

Then 



□ 



1 /•(0+) 



1 1 



■Kt {2Viy 



(66) 
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Proof. We compute initially the above integral in the particular case m — for 
which by the reparametrization z = w'^ and an elementary transform we obtain: 

r'(0+) _ -| pc+ioo 

' c—%oo 



1 r^^^> r- 1 r 

— \ (v^)-i e-«^ e^* = — 2 / 



27ri 

It is easy to see that the above relationship is differentiable any number of times 
on a under the sign of the integral, hence 

\ ^) I I rZ\m-\ „-oiJ^ „zt J ^ 



2m 



e-«^e^*d0 = ^^e-" Z^^*' = (67) 



1 5'".,-(a/(2Vt))= 1 1 ^" 



aa™ (2^^)"^ du" 

_J_j-ir 

^^rt(2^/^)- '"Ux/i/ 
where the last equation has been written by virtue of the well-known relationship 

(67) implies 

' x/^(2\/i)'» \1^t) 

and this completes the proof. □ 

Corollary 2. 

27r^i_o„ ^ V2^(2V2t)™ \2y/2t J 

According to Batcman and Erdelyi (1954) 5.6 (12) we can write 
1 /■(°+' e""^ e""'/^''*) 



27rz 



-OO + V^t '\2^t ^ , 



where a, /? are real parameters. 
Introduce the notation 



$o(a,/3,t) := - /?e«/'+/'^* er/( ^ + /3Vi 

$™(a,/3,i) := (-l)"|n.cl>o(a,;3,f), 



(68) 



(69) 



where m = 1, 2, 

Remark. Observe that repeated differentiation on a of the function ^o{a,p,t) 

results again in the error function of the same argument, er/^^^ + (3\/i^ , which 

is multiplied by a certain exponential function with some exponential functions 
added. This follows from the definition 

2 

erf{x) = — I e""' du (70) 
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and elementary rules of differentiation. The functions ^m{o., P,t) are computable 
recursively by means of symbolic differentiation; for computations they are as good 
as elementary functions. 

Proposition 8.4. The integral in (70) as a function of the parameter a is differ- 
entiable under the sign of the integral any number of times, hence 



where m = 1, 2, 



Proof. The last integral in (71) absolutely converges for any m. The proof of (71) 
follows from (70) by induction on m and by a standard application of bounded 
convergence theorem. □ 

Corollary 3. 

1 /•(°+) (^/i)"a-v^ 



27rj 



/ ^V^j a e^'Uz = ^m{ina,k,2t). 

J-oo Vz + k 



The following theorem completes Laplace inversion; it is the main theorem of 
the present article. 

Theorem 11. The density distribution f{X,t) of functional (1) is determined by 
the following absolutely convergent series: 

^ ' 2XV2X ^n! ■'"^ ^(2A)"' ^ ^ 



where 
and for n = 1,2,... 

fn{X,t)= (74) 



where a = 8Ae Hm{x) are Hermite polynomials, ^m{c(, P,t) are defined by (69), 

(n 



the coefficients dm are defined by (54) - (55), and the coefficients a^"'' are defined 



in Theorem 7. 

Proof. The expansions above and their absolute convergence follow from Theorem 
6 and Proposition 8.2 by substitutions indicated in Corollary 2 and Corollary 3. 

□ 

Remark. Absolute convergence of these series can be also seen from the estimates 
of their remainders given in the next section, see Theorem 12 and Corollary 4. 
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9. Asymptotic Error Bounds 



Estimates similar to those used in Section 7 readily lead to accurate error bounds 
that allow computations by series given in Theorem 11 with guaranteed accuracy. 
This matter, however, exceeds the limits of the present article. In practice the 
required accuracy is usually achieved by other means of numerical analysis. For 
these reasons we content here with asymptotic estimates of error bounds. 
The following theorem estimates the remainders of the scries in (72) (74). 
Let i?Ar(A,i), QM(A,i) and Q^J^\X,t), where k = 1, . . . , [(n — l)/2] + 1 be these 
remainders: 



°° I 1 



n=N 

oo 



QM{X,t) = 



dn 



Hm 



{2\Y 
/ Ina \ 



m=M 



Theorem 12. For fixed A, t > and N,M ^ 00 holds: 

(i) 



\RN{\,t)\=0 



(ii) 
(in) 



(2A)^iV!y' 
|Qm(A,OI =0(M-(i-^)^/2), 
|Q(^)(A,t)|=0(M-(i-W2), 



where k = 1,. . . , [(n — l)/2] + 1 is fixed 
Proof. Proof (i). 



\RN{\t)\ 



°° 1 1 °° 1 1 

E.;^^"(^'^)(2aF <E,;^i^"(^'^)im 



n=N 



' (2A)^ 



where by (36), (37) and according to (29) 



l/n(A,t)| 



27ri A. 



c+ioo 



w„{\,q) e«* dq 



1 



< 



C— ICX) 

c+ioo 



27ri Jc—ioo 



-v/2^/2 r(^+n) 



r(^/2^ + n+l) 



e«* dq 



< 



|(2A) 



-V2?/2 



2tt 

:= C{\,t) < 00, 

where C(A, t) does not depend on n 
This implies 



|2V^- 



^1 



r(y2^/2 + i/2) 



\dq\ :-- 



\R.iX,t)\ < C(A,t) £ 1^ = ^^0(1) = O(^). 
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Proof (ii). By (63) and estimates (56) and (61) we have 

oo „(0+) 

|Qm(A,OI= E I^-I ^ / {^zT-^l^ a-^^ e^^Uz 



0{M- 



-(l-e)M/2 



Proof (iii). In a similar fashion to Proof (ii) by (63) and estimates (56) and (61) 
we have 



\Q'-M\\t)\ 



'^■^ 



oo 
m=0 



2m 



1 



m 



-oo + ^ 

(0+) 



< 



-sfz 2zt 



2-111 ./-^ \/z + k 



dz 



□ 



Corollary 4. t/ie series in Theorem 11 converge absolutely. 



These error bounds show that for A's not too smah the convergence rate of ah the 
series in question is exceptionally good; for the small A's, however, an asymptotic 
expansion is available to fill this gap. This expansion will be discussed in a paper 
to follow. 

10. Computations 

The above series and error bounds provide effective means for computations 

of f{X,t), as illustrated on the example of the density distribution /(A, 1) of the 
functional e^^^^^^ds, see Table 1. A graph of /(A, 1) and further information on 
this distribution are available via www.tolmatz.net. 
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14 
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1 


15 
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1 
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16 
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1 
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1 


17 


0.2949056229 


1 
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0.18 





3696381349 


0.68 
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18 
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1 
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Table 1. The density function /(A, 1). 



